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ABSTRACT 



We present a multi-wavelength study of a 3.6/im-selected galaxy sample in 
the Extended Groth strip. The sample is complete for galaxies with stellar mass 
>lO9-^M0and redshift OA < z < 1.2. In this redshift range, the IRAC 3.6//m 
band measures the rest-frame near-infrared band, permitting nearly unbiased se- 
lection with respect to both quiescent and star-forming galaxies. The numerous 
spectroscopic redshifts available in the EGS are used to train an Artificial Neural 
Network to estimate photometric redshifts. The distribution of photometric red- 
shift errors is Gaussian with standard deviation ~0. 025(1 + z), and the fraction 
of redshift failures {>3a errors) is about 3.5%. A new method of validation based 
on pair statistics confirms the estimate of standard deviation even for galaxies 
lacking spectroscopic redshifts. Basic galaxy properties measured include rest- 
frame U — B colors, B- and K-band absolute magnitudes, and stellar masses. 
We divide the sample into quiescent and star-forming galaxies according to their 
rest-frame U — B colors and 24 to 3.6 //m flux density ratios and derive rest K- 
band luminosity functions and stellar mass functions for quiescent, star forming, 
and all galaxies. The results show that massive, quiescent galaxies were in place 
by z ~ 1, but lower mass galaxies generally ceased their star formation at later 
epochs. 

Subject headings: cosmology: observations — galaxies: redshift — galaxies: Survey — 
galaxies: mid-infrared 
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INTRODUCTION 



The central puzzle of galaxy evolution is how baryons fall into galaxies and turn into 
stars. Understanding the mass assembly history of galaxies is a critical element for solving 
this puzzle. Galaxies can be generally divided into two populations: those actively forming 
stars and those that are quiescent with relatively little star formation. Observationally, 
these two types of galaxy have differing colors corresponding to different regions in 



color-magnitude diagrams (iStrateva et al. 



2001 



Hogg et al. 



20031 ). The star- forming 



galaxies are said to occupy the "blue cloud" (BC) while the passive galaxies occupy the "red 
sequence" (RS). Similar bimodal distributions are also seen in their naorphologies, spec tral 



types, metallicities, and star formation rates (IMadgwick et al. 



2003 



Kauffmann et al. 



2003al Jbl). Morphologically most galaxies in the RS are elliptical/SO galaxies, and galaxies 



in the BC are disk or irregular /pe culiar galaxies. This bimodality is observed over a very 



wide redshift range u p to z < 2.5 (IBell et al. 



2004 



Willmer et al. 



2006 



Faber et al. 



2007 



Brammer et al. 



20091 ). In a color-stellar mass diagram, red galaxies appear to be much 



more massive than blue ones. 



Formation of massive red galaxies cannot be accounted for 



galaxy evolution 


path. Many aut 


lors 


Toomre & Toomre 


1972 


Kauffmann et al. 


1993; 


Cole et al. 


1996; 


Hopkins et al. 


200 


5) su 



by any simple blue-to-r ed 



Toth fc Ostriker 



1992 



20061 ) suggest that galaxy merging 



is the only way of forming massive red galaxies. There are probably more than two merging 
scenarios: two blue galaxies can merge to trigger a burst to star formation and eventually 
form a massive red galaxy, or two blue galaxies can evolve to red galaxies, then undergo 
a dry merging to form a massive galaxy. A census of red and blue galaxies at different 
redshifts is needed to constrain how galaxies migrate from the BC to the RS. 



Redshifts are essential 



surveys fiLilly et al. 



1995 



for studying galaxy populations 



Cowie et al 



1996 



Lin et al. 



. Early spectroscop ic galaxy 



1999 



Cohen et al. 



2OO2I) usually 
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obtained several hundreds of redshifts. These surveys concluded that strong galaxy 
evolution seen in faint S-band galaxy counts was mainly d ue to blue galax y evolution 
since z = 1 with red galaxies showing little or no evolution. 



Imetal 



(120021) constructed 



luminosity functions for elliptical galaxies in 0.5 < z < 1 with the DEEPl redshift sample 
and argued that they saw only passive luminosity evolution but no change in number 
density for elliptical galaxies since z = 1. These early redshift surveys are too small to have 
significant statist ics for massive red galaxies, which can also suffer a large cosmic variance 



(ISomerville 



20051 ). Recently several la rge optical redshift s urveys have b een carried out o n 



8 m class tele scopes including VVDS ( 



and DEEP2 (IFaber et al. 



2007 



Le Fevre et al. 



Davis et al. 



20051), zCOSMOS dhillv et al.ll2007h . 



20071 1. each with more than 10000 redshifts 



available for galaxies at < 2; < 1.5. These redshift surveys are highly complete for blue 
galaxies up to 2; = 1.4 because the [O II] A3727 line can be observed. Measuring redshifts 
for red galaxies is more difficult because there are only absorption lines available in most of 
their spectra, and at high redshifts, only luminous red galaxies can have redshifts measured 
in these existing redshift surveys. Therefore optical redshift samples suffer severe bias 
against red galaxies at high redshifts. 

Most studies of evolu tion of galaxy populations have u sed optically-selected 



samples. Several groups (jBell et al. 



2006 



Faber et al 



20071 ) used photometric or 



photometric/spectroscopic combined redshift samples to study evolution of red galaxies 
and found that the number density of red galaxies increases from z = 1 to 2; = 0. This is 
consistent with a migration of galaxy population from BC to RS. However, as mentioned 
above, optically-selected samples are biased against red galaxies at 2; ~ 1, where even I 
band probes the rest-frame UV. 



Near-IR selected galaxy samples permi t detection of red ga 



or at higher redshifts (IBrammer et al 



2009 



van Dokkum et al 



axies with lower luminosities 



20091 ) than optically-selected 
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samples. The rest-frame -ft^-b and light of a 



traces underlying stellar mass( ICowie et al. 



galaxy comes mainly from M and K stars and 



1996 



Huang et al. 



1997 



Bell fc de Jong 



200l|). 



Therefore a rest-frame i^'-band selected galaxy sample is equivalent to a mass-selected 
sample unbiased against either red or blue galaxies, permitting fair statistics for both galaxy 
populations. However, early ground-based near infrared (NIR) surveys have been ineffective 
because of small format IR array camera and high NIR sky background. Large-format 



coverage and faint 



2003 



Bundy et al. 



2006 



Tavlor et al. 



2009 



f» 23-25 mag ( 


Cima 


;ti et al. 


2002 


Franx et al. 


van Dokkum et al. 


2009 




McCracken et al. 





20101 ). Most deep i^-band sur veys focus on forrnation and evolution of passive red galaxi es 



with > 10^^ Mr;, at ^ ~ 2 ( 



van Dokkum et al. 



2009 



~ 2 ( 


Franx et al. 


2003; 


McCracken et al. 


2010) 



Taylor et al. 



2009; 



Brammer et al. 



2009 



20101 ). On the other hand, loc al i^-band galaxy 



surveys detect e arly type galaxie s with stellar r nasses as low as 10 M© ( iCole et al 



Bell et al. 



2003f ). More recently. 



Baldry et al 



2001 



( 120121 ) estimated the local stellar mass 



function down to M^,=1O''M0, and found that the local stellar mass f unction has a yery 



Bundy et al. 



(120061 ) 



steep slope at the lower mass end. Most relevant to the present study, 
carried out i^"-band imaging in the Extended Groth Strip region with Kab < 23. The 
limiting stellar mass for their fT-selected sample is ~1O^°'^M0 at z ~ 1. Formation and 
evolution of lower mass red galaxies since z = 1 has not yet been well studied. 

The InfraRed Array Camera (IRAC) on board the Spitzer Space Telescope can perform 
very fast scanning of the IR sky with great depth. The IRAC 3.6 and 4.5 fim bands probe 
the rest NIR for galaxies at z ~ 1. There are several large, deep IRAC surveys. Samples 
obtained in these surveys are usually confusion limited at /jy ~ 1 /iJy. The limiting flux 
of 1/iJy ([3.6] = 23.£lj) corresponds to an absolute i^-band magnitude of Mk = —19 for 



4n this paper, the notation [w] means the AB magnitude at wavelength w in /xm. Some 
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a galaxy at ^ = 1 . A galaxy with [3.6] = 23.9 mag has a typical i?-band magnitude of 26 



( iBarro et al 



201ll ). too faint for spectroscopy even with 8m class telescopes. Photometric 



redshifts are essential in studying faint mid-infrared (MIR) selected galaxies. 

Using photometric redshifts is inevitable in studying red galaxy evolution. With large 
format CCD cameras available on 6-8 m class telescopes, deep multi- wavelength photometric 
data become available. These data permit measurement of accurate photometric redshifts 
for much larger samples of faint galaxies than can be obtained spectroscopically. A difficulty 
is that for galaxies at higher redshifts, the dominant features in galaxy SEDs, for example 



the Balmer/4000A break, shift out of the visible b and. Including intermediate- 



3and 



2003 



Ilbert et al. 



20091) 



photometry can improve the photometric redshifts (IWolf et al. 
and accuracy can be further improved by extending the wavelength coverage to the NIR. 
Recently the NEWFIRM intermediate band survey, including 3 bands in the 1.2 /xm window 
and 2 bands in the 1.6 /im win dow, was completed w ith the NO AO 4 m telescope in search 



of passive galaxies at z > 1.5 f lWhitaker et al. 



201ll ). The relative redshift uncertainty 



(Jz/il + z) for their photometric redshifts is 1-2% for the K-selected galaxies. 

To extend SED fitting to the Spitzer IRAC bands, photometric redshift estimation 
must overcome a new challenge: variation of point spread function ( PSF) from ultraviole t 



to the mid-infrared. Traditional 



Brodwin et al 



2006 



Ilbert et al. 



y, photomet ric redshift estimation (iRudnick et al. 



2006 



2001 



20091 ) is performed by fitting a set of templates 



to an observed galaxy spectral energy distribution (SED) obtained in a multi-wavelength 
photometric survey . PSF variation in different bands may distort the observed SEDs 

20051) ■ and anyway we have poor knowledge of galaxy SED templates in the 



(iLabbe et al 



MIR ( 



Lu et al. 



2003 



Brodwin et al. 



20061 ). Artificial neural networks (ANNz, CoUister & 



Lahav 2004) is an alternative redshift estimation method. ANNz has several advantages for 



other papers have used this notation to mean Vega magnitudes. 
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MIR-selected galaxy samples: it does not require any knowledge of galaxy SED templates, 
and the resulting photometric redshifts are not affected by photometry offsets caused by 
PSF variation or any other systematic offset in different bands. ANNz uses a calibration 
sample to set up an empirical relation between spectroscopic redshifts and observed galaxy 
properties. Input parameters can be anything including color, flux, size, and even digital 
morphological parameters. All this method requires is that input and calibration samples 
have the same filter set for photometry and cover same input parameters and redshift 
ranges. Though ANN^ canno t measure any redshifts beyond the range of the calibration 



sample. 



Brodwin et al. 



( 120061 ) showed that ANN2; photometric redshifts for a MIR-selected 



sample have much smaller scatter than those derived with the SED template methods. 

This paper presents a mult i- wavelength study of a complete 3.6 /xm-selected galaxy 
sample aX z < 1.2 in the Extended Groth Strip region. Photometric redshifts for this sample 
are derived using ANNz. §2 describes the sample selection, and §3 explains the photometric 
redshift derivation with ANN2 including the uncertainties. §4 gives results for the blue and 
red galaxy populations at 0.4 < 2; < 1.2 and their properties including the rest K-band 
luminosity function and color-magnitude relation. §5 summarizes results. Source distances 
are based on standard flat ACDM cosmology with Hq = 70 km s~^ Mpc and Qm = 0.30. 



2. The 3.6 /xm Selected Sample in EGS 



The All-wavelength Extended Groth Strip International S urvey (AEG 



20071) is a large mu lti-wavelength survey in bands from X-ray (iLaird et al. 



( Ilvison et al 



Davis et al. 



20091 ) to 20 cm 



20071 ). The original Chandra survey covers the whole strip (2°x 15') with 



200ks per pointing, and later the central 45'xl5'strip with a deeper exposure of 800ks per 
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pointing. The x-ray limiting flux densities □ for the 200ks imaging in the whole strip are 
5.3xl0^^^erg cm ~^ s~^ in the soft ( 0.5-2 keV) band an d 3.8xl0^^^erg cm~^ s^^ in the hard 



(2-10 keV) band flLaird et al 



20091). 



Ivison et al 



( 120071 ) performed a deep radio imaging at 



1.4GHz in the EGS region with the 5a limiting flux density of 50/iJy. Both X-ray and radio 
can help to identify AGNs. More recent data include Herschel Space Telescope FIR/subrn m 



imag i ng in the EGS as 



2012: 


Smith et a . 


2012) 


mJy( 


Oliver et al. 


2012). 



part of the HerMES survey (Oliver et al. 



201G 



Roseboom et al. 



20121) . The SPIRE 5a limiting flux densities are 13.8, 11.3, and 16.4 



20121 ). The angular resolution for the SPIRE images are rather poor, 18", 



25", and 36"at 250, 350, and 500/im, sources a re blended in the EGS SPIRE images. Deep 



20071 ) with the 5cr limiting 



24 /im imaging is available in the EGS region (iDickinson et al. 
flux density of = 20 /iJy, providing a good prior to identify the blended Herschel SPIRE 
sources. The actual FIR photometry catalog in the ACS covered a rea was extracted using 



the MIPS 24/xm catalog as the prior with the method described in 



Rigopoulou et al. 



Roseboom et al 



fl2012h 



(120 12l ) derived FIR luminosities for Herschel sources in the region by 



fitting templates to their IR SEDs. 



2003; 


Faber et al. 


2007 


Cooper et al. 


2011, 


2012) 



E 



The AEGIS datasets include one of the 



deepest IRAC surveys, covering 10' x 2° with >9ks exposure time. 



data reduction, and photometry were 



more aggressive source extraction by 



presented b y 



Zheng et al. 



Barmbv et al. 



'he 



RAG observations. 



fl2008[ ) and later with 



mm . The latter authors found 80% 



completeness at limiting flux densities of 1.8 and 2.0 /xJy respecti vely at 3.6 um a nd 4.5 /im 
compared to ~60% completeness at these flux densities found by Isarmby et al.l . For the 



^The lir aiting flux density is defined as the flux to which at least 1% of the survey area 
is sensitive (ILaird et al.ll2009l ). 

^http://deep. berkeley.edu/DR3/ 
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current paper, w e have restricted our sample to the HST ACS imaging region (10'xl°, 



LotzetaL 



20081 ). where there are generally deeper visible and NIR photometric data 



available. Our prel iminary sana ple includes all objects in this region with flux density 



/s.e > 2 fiJy in the 



Zheng et al. 



catalog. This limit selects galaxies over a wide redshift 



range and also a few stars. There are 23327 objects in the preliminary sample. 

The visible and NIR photometric data for the sample are summarized in Table 1. 
There are 18 bands from m to 8 fim that are useful for photometric redshifts. Figure [1] shows 
a histogram of the number of photometric points available per object. About 96% of the 
objects have photometry in at least 5 bands, and 62% of the objects have at least 15-band 
data. As Table 1 shows, however, some of the visible and NIR bands are redundant. For 
example, there are four versions of i?-band imaging with various depths. In redundant 
bands, we used a cx-weighted mean to combine photometry to give a single photometric 
value with maximum signal-to-noise ratio in each band. After combining redundant bands, 
there are 12 wavelengths available in the u to 8 M m range. In the y isible range ugriz, the 



Canada France Hawaii Telescope Legacy Survey flllbert et al. 



2006 



, CFHTLS) provides a 



well-calibrated photometry system, and we transformed all photometry in this wavelength 
range into the CFTHLS system. 

This study focuses on galaxy population at 0.4 < z < 1.2 to take advantage of the 
large number of spectroscopic redshifts available. Therefore the preliminary all-z sample 
needs to be refined by rejecting foreground stars and background galaxies at z > 1.2 . A 



set of simple color crite ria can separate the background galaxies. 



Huang et al. 



(120041 ) and 



Sorba fc Sawickil (120101 ) show that the IRAC [3.6] — [4.5] color is a function of redshift 



(Figure [2]): galaxies with [3.6] — [4.5] < are at z < 1.25, and galaxies with [3.6] — [4.5] > 
1.250 



are at ;z > 



This is because 1.6 /xm bump in galaxy SEDs shifts through the IRAC 



"^The exact redshift for galaxies with [3.6] — [4.5] = is weakly dependent on galaxy 
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bands at z > 1. Several groups have used this SED feature to select galax ies at ^: ^ 2 for 
Spitzer IRS spectroscopy with a nearly 100% successful identification rate (IWeedman et al. 



2006 



Farrah et al. 



2007 



Huang et al. 



2009 



Desai et al 



2009 



Fadda et al. 



2010h . Some 



rest-frame UV-selected galaxies at 2 < 2 < 3 such a s BM/BX ga. 



axies and LBGs, howeve r, 



2005 



Rigopoulou et al 



are t oo young to have developed the 1.6 /im bump ( iHuang et al. 
20061 ). and these galaxies may have [3.6] — [4.5] < even with 2; > 1.2. We therefore need an 
additional criterion to reject these young blue galaxies. We used an optical-MIR color-color 
diagram {R — I vs R— [3.6], Figure |3]) for classification. Figure [3] shows that stars, galaxies, 
at 2; < 1.25, and galaxies at z > 1.25 occupy different regions. Thus the three criteria for 
selection of the 3.6 fim galaxy sample are: 



[3.6] - [4.6] < 

R-I > 0.2{R- [3.6]) + 0.015 
R-I < 0.4(i?- [3.6]) + 0.5 



(1) 

(2) 
(3) 



Criterion 1 removes background galaxies with a normal 1.6 /im bump, criterion 2 rejects 
young, blue galaxies at z > 1.25, and c riterio n 3 rejects stars. Criterion 1 is not as stringent 
as that proposed by 



Sorba fc Sawickil (120101 ) for the same purpose, thus may reduce the 
bias by the color cut, but will include galaxies at z>1.3 in the preliminary sample. These 
galaxies will be rejected from the sample after we apply Criterion 2. Criterion 2 is very 



similar to t 



l e BzK selection for galaxies at z~2, but used here to reject these galaxies. 



Ilbert et al. 



( I2OO9I ) used their SED-fiting photometric redshift sample to demonstrate that 



galaxies can be well separated at z=1.3 in the same R— 1 vs R— [3.6] diagram. 



spectral types w i th re d galaxies having z = 1.25 and blue galaxies having z = 1.50. 
Sorba fc Sawickil (120101 ) proposed a more stringent color-dependent criterion to separate 
galaxies at z=1.3 : [3.6]-[4.5]=0.12([3.6]-[8.0])-0.07 
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By combining both Criterion 1 and 2, we can effectively reject galaxies at z>1.25 
while reduce the bias caused by applying only one color criterion. The final sample for 
photometric redshift estimation includes objects in the HST/ACS area of the EGS meeting 
the above three criteria and with [3.6] < 23.15. Figure H] confirms that our sample selection 
does not have any strong bias against either quiescent or star forming galaxies. There are 
12890 galaxies in the final sample. 

Figure [5] illustrates successive application of the color criteria. In our preliminary 
sample, 95% of the objects have both R and / band detections and can be selected by all 
three criteria. Of 1297 sources missing either R or /, 1185 have 4.5 /im counterparts which 
we will be able to classify with relation 1: 387 with [3.6] — [4.5] <0 are included in the sample 
and 798 with [3.6] — [4.5]>0 are reject. The remaining 112 sources cannot be definitively 
excluded and are kept in the sample, but their photometric redshifts will be more uncertain 
than sources having better wavelength coverage. 

Our selection does not particularly discriminate against AGNs, but AGNs with 
power-law SEDs in the mid-infrared band. Objects with strong hot dust emis sion are calle d 



' pow er-law" objects with very much red color of [3.6] — [4.5]>0 even at z<l (jStern et al. 



20051 ). Thus their 3.6/im flux densities for these power-law objects a re mainly from hot dust 



Huang et al. 



( 120071 ) found that number of 



emission, but no longer trace their stellar mass, 
power-law objects is less than 1% of >L* galaxies in an 8/im selected sample. The 3.6/im 
band is less sensitive to power-law objects than the 8yum band, thus their fraction is even 
lower in the 3.6/im selected sample. These objects will be rejected from the sample by 
Criterion 1. 
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Photometric Redshift Estimation 



The ANNz method uses a subsample with known redshifts as a training set to build 



up an empirical relation between measurable parameters and redshifts a nd then app 



relation to the rest of the sample to derive their photometric redshifts (IFirth et al. 



Brodwin et al. 



ies this 



2003 



20061 ). Almost any measurable parameter of galaxies including magnitude, 
surface brightness, color, size, and morphology can be used as an input node for ANNz, 
but magnitudes at a range of wavelengths are the most common. ANN2; requires that the 
training set should cover the same ranges as the photometric sample in all input parameters. 
In other words, ANNz cannot derive photometric redshift for a galaxy with input parameter 
values beyond the range of its training set. The spectroscopic redshift sample has a limiting 
magnitude of i? < 24.1, but a significant number of objects in our sample are fainter than 
this limit. Using magnitudes as input parameters to train the ANN2; will allow it to derive 
photometric redshifts only for galaxies with R < 24.1. We therefore used colors with 3.6 fim 
photometry as the color base as the input nodes in the ANNz estimation. There are 11 
non-redundant colors for each galaxy as shown in Table 1. Faint galaxies {R > 24.1) at 
z < 1.25 exhibit the same color range as the galaxies with spectroscopic redshifts shown 
in Figure El and therefore the program will be able to drive photometric redshifts for the 
whole sample. 



ANNz network architecture can be described as A^input : : A'2 : : Nm '■ A'output, 

where A^input is the number of input colors for each object. For the present data, A'input = H, 
the number of non- redundant colors. A^i, N2, .... and are the number of nodes in each 



ANNz hidden layer for a total of m layer s, and Armtp n t, is t 



photometric redshift, usually set to be 1. 



Firth et al. 



le number of sets for output 



(120031 ) did extensive experimenting 



with ANN2; to determine the optimal number of hidden layers for the best photometric 
redshifts and found that a multi hidden layer architecture yields photometric redshifts 
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with a lower (y{zp — Zg) than a single hidden layer architectu re. A large numbe r of layers, 



Firth et al. 



mom adopted 



however, do not significantly improve photometric redshifts. 
6:6:6:6:1 as their fiducial architecture. We ran the ANNz with various architectures for our 
sample and arrived at the same conclusion. The final redshifts use 11:8:8:8:1 architecture, 
meaning 11 input colors, 3 hidden layers with 8 nodes in each layer, and 1 output. 

The accuracy of the AN Nz photomet r ic red shifts is partly determined by the number 



of objects in the training set. 



Firth et al.l ( 120031 ) used a sample of 20000 objects to show 



that a larger training set can yield more accurate photometric redshifts, but increasing 
training set size beyond retrain = 1000 did not improve the accuracy significantly. There 
are ~5400 spectroscopic redshifts available in the sample. We divided the spectroscopic 
redshift subsample into 3 data sets: a training set, a validation set, and an independent 
validation set. The first two sets are required by ANN2;, and the last one is used to give an 
independent photometric redshift validation. 

Any bad entries, either bad photometry or bad spectroscopic redshifts, should not 
be included in the training set. Including these outliers will yield large scatters for the 
resulting photometric redshifts. We ran the ANN2; estimation twice. After the first run, 
we compared photometric redshifts with spectroscopic redshifts in the sample and fiagged 
those with \zp — Zs\/{1 + Zg) > 3(7 For a second training run, we updated the training and 
validation sets omitting fiagged objects, about 3.5% of the initial sample. This iteration 
improved the accuracy of the photometric redshifts significantly as shown in Figure [7] and 
Figure El 

It is critical to assess the accuracy of photometric redshifts before using them. 



^We define dp, the calculated uncertainty of the photometric redshifts, as the dispersion 



of SZp = {Zp - Zs)/ (1 + Zs). 
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Traditionally objects with spectroscopic redshifts are used to evaluate photometric redshift 
errors. ANNz requires one validation set of objects, and we created an independent 
subsample with spectroscopic redshifts, none of which is used for the ANN2; training, for 
an additional validation. The comparison between photometric and spectroscopic redshifts 
in Figure [7] permits a direct estimation of photometric redshofts errors. The photometric 
redshifts derived in the first run show much more systematic deviations caused by the bad 
redshifts in the training set. Figure [8] shows a Gaussian distribution of redshift errors, 
implying a well-defined random uncertainty characterizes most of the photometric redshifts. 
Both the ANN2; and independent validation sets yield similar dp ~ 0.025 in the second 
run, again with about 3.5% outliers. The training set cxp cannot be used for an uncertainty 
estimate because it was deliberately minimized by the calculation. The outlier galaxies have 
a median 3.6 fim flux density of 12 /iJy, much higher than the 6.5 /iJymedian flux density for 
the whole sample. The higher average 3.6 fim flux densities outlier galaxies may have IRAC 
fluxes contaminated by nearby objects, or two sources are unresolved by IRAC but resolved 
at shorter wavelengths. The IRAC 3.6 /im image has a 2" angular resolution, whereas the 
CFHT i?-band images from which the DEEP2+3 sample was selected have resolution of 
~1". Contamination or source blending are therefore expected to occur. 

Galaxies with spectroscopic redshifts tend to be the brightest galaxies in any sample, 
and this is especially true in our 3.6 yum-selected sample. Galaxies that are faint in visible 
light have relatively larger photometric errors and therefore may have larger errors in 
their photometric redshifts. Validation with the spectroscopic redshift subsample may not 
characterize the photometric redshift scatter for the fainter galaxies, and it is essential to 
assess photometric redshift accuracy for these faint objects. 

We propose a new method to estimate the photometric redshift standard deviation 
for the whole sample. The method is based on galaxy-galaxy pair statistics. A substantial 



- 16 - 



fraction of galaxy visual pairs are actually physically related, and the redshift offset for 
a real pair should be near zero. However, the calculated photometric redshift offset for 
a real pair may be nonzero if either redshift is erroneous. We define a visual pair as 
any two galaxies within an angular separation of 30" and define the redshift offset as 
6zss = {zsi — Zs2)/{1 + Zsi), where Zsi and Zs2 are spectroscopic redshifts for two galaxies of 
the pair. The 30" angular distance is corresponding to a distan ce of 160-240 kpc for galaxies 



Linetal 



(120071) to study 



LinetaL 



tooi ) 



at z=0.4-1.2, similar to the criterion of d<300kpc adopted by 
galaxy pairs and mergers in the same redshift range in the same field, 
found that these galaxy pairs have a very wide range of star formation activity. Galaxy 
pairs with shorter separation (~50kpc) likely have higher star formation activity, otherwise 
they do not show any bias towards any type. 

Figure M show the histogram of 6zss for all visual pairs with both galaxies having 
redshifts in the DEEP2+3 sample. The histogram shows two components: a spike at 
6zss = for real pairs and a broad distribution for projected pairs that are physically 
unrelated. The expected 6zss for projected pairs can be easily generated using Monte-Carlo 
simulation and subtracted from the distribution for all pairs. The residual spike is the 
Szss distribution for real pairs. As expected, both of these velocity components are small 
compared to the expected uncertainty in the photometric redshifts. 

Having a large number of spectroscopic redshifts available allows pair statistics to 
give a good estimate of the uncertainty in the ANN2; redshifts. The pair sample tested 
consists of galaxies separated by <30" where one galaxy has a spectroscopic redshift 
and the other only a photometric redshift. The normalized redshift offset is defined as 
^^sp = (zp — Zs)/{1 + Zg), where Zp is the photometric redshift and Zs is the spectroscopic 
redshift of the respective galaxies. The histogram for all pairs (Figure [TOl) shows two 
components, as with spectroscopic pairs. Subtracting the simulated distribution for 
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unrelated pairs from the real pair distribution gives the 5zsp distribution for physical pairs. 
This distribution shows a Gaussian shape (Figure [TU]) with (Xsp = 0.024. This is very close 
to (Tp, the standard deviation of the photometric redshifts estimated from the difference 
between spectroscopic and photometric redshifts for the same galaxies. 

The pair method can be used directly to validate the photometric redshifts for faint 
galaxies with R > 24.1. Figure [TT] shows the results. 6zsp for these galaxies also shows a 
Gaussian shape with asp = 0.024. The pair statistics thus confirm the validity of the ANN^ 
photometric redshifts even for galaxies fainter than the limiting magnitude of the training 
set. 

Overall, dp of the photometric redshifts for the whole 3.6 fim sample is very close to that 
of recent state-of-art SED-fitting photometric redshifts for optically-selected s amples, whi ch 



often have more 


lan ^ 


2003 




Ilbert et al. 


2006 



lan 20 bands of intermediate and broadband photometry (jWolf et al 



2009 



Whitaker et al. 



201 ll ). Comparing with those photometric 



redshift estimations, the ANN2; redshifts have well-defined and consistent error distributions 
for both optically bright and faint galaxies but are inherently limited to z < 1.25, the 
maximum redshift of the training set. 



4. Galcixy Populations in the 3.6 /xm Selected Sample 



4.1. Color Classification 



With redshifts known, deriving rest frame colors and magnitudes for sample galaxies 
is straightforward. All types of galaxies have very similar SEDs at NIR wavelengths. The 
i^-band absolute magnitudes fo r galaxies can therefore be calculated with the 3.6 /im flux 



densities and K-correction from 



De Propris et al. 



(I2OO7I ). The absolute magnitude range for 



the sample is —19 < Mk < —25. Rest frame U — B colors for the sample were derived from 
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2006 



Faber et al. 



V — I col ors via a set of lo cal galaxy SED templates (iKinney et al.lll996l — IWillmer et al 
and 



20071 took a similar approach) . Rest frame absolute B magnitudes 



Mb were derived in the same way. The RS and BC (Ivan Dokkum et al 



with the 



Willmer et al. 



2000l ) were divided 



( I2OO6I ) criterion: 



{U - S)o = -0.032(Mb + 21.52) + 0.454 - 0.25 + 0.832 



(4) 



using the same U and B magnitude system they used but adding 0.832 to convert from the 
Vega to the AB system. Figure WI\ shows a series of color-magnitude (CM) diagrams for the 
3.6 /xm-selected sample. 

Two types of galaxies have colors placing them in the RS: dusty sta r-forming galaxies 



and quiescent galaxies with little or no star formation (ILabbe et al 



20051 ). Current studies 
use e ither a rest-frame visible color-color diagram or IR emission to separate these types 



(e.g. 



Brammer et al. 



2009 



Williams et al. 



2009; 



Wang et al. 



2OI2I). 



Williams et al. 



mm 



showed that the types can be separated in a f/ — V ys V — J diagr a m, w. 



have redder V — J colors than quiescent galaxies. iBrammer et al 



lere dusty galaxie s 



mm and 



Wang et al 



( 120121 ) used 24 /im emission to identify dusty star forming galaxies. The 24-to-3.6 /xm flux 
density ratio, f2il is closely related to the specific star formation rate and should 
therefore provide a good way of distinguishing passive from dusty star-forming galaxies. 
Figure [13] shows that f^i/fss is roughly constant for 24 /im sources with BC U — B colors, 
but this ratio decreases to below 1 for 24 /im sources with RS colors. Nearly all 24 /xm 
galaxies in the BC have f2i/ fz.e > 1, whereas more than 50% of 24 /xm galaxies in the 
RS have f24/ fss < 1- All i/ersc/ie//SPlRE sources in the RS also have f24/ fss > 1- 
These sources must be star forming to be detected with SPIRE. We have therefore used 
f2i/ fz.Q = 1 to separate dusty from passive galaxies in the RS: galaxies with /24//3.6 < 1 
are taken to be passive galaxies with 24 /xm residual emission, and galaxies with f2i/ fz.& > 1 
in the RS are taken to be dusty and star forming. This criterion works because most RS 
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galaxies in the sample are very luminous with /a.g > 20 /iJy and therefore RS galaxies not 
detected at 24 /xm must have f2i/ fzs < 1- 



Galaxies with colors near the dividing line between the RS and the B C are called 



valley" objects. Most green valley galaxies are also MIPS 24 yum sources (ICowie et al. 



Kartaltepe et al. 



green 



2008 



20101) • An IR-selected sample has many more green valley objects than a 
visible-selected sample, and this is the case for our 3.6 /im-selected sample. Figure [12] shows 
many MIPS 24 yum sources in the green valley, and some are well into the RS color zone. 
On the other hand, most galaxies with no MIPS 24/im detection have rest U — B colors 
well away from the green valley, and the rest U — B colors should accurately classify these 



galaxies as passive or star-forming. 



4.2. Star Formation and Stellar Mass Assembly in the 3.6 /im Sample 



Stellar mass for the 3. 6 /xm sample can be deriy ed from stellar population models 



In practice we used BC03 f iBruzual fc Chariot 



20031 ) models with solar metallicity and 



Salpeter initial mass function (IMF). The model template set was generated with a grid 
of stellar masses, E{B — V), ages, and e-folding times for star formation history. Stellar 
masses for the sample galaxies were determined by fitting their observed SEDs to the 



mode 



templates. 



IMF flBundv et al 



(IBell et al. 



he derived stellar masses can be converted to those with the Chabrier 



2006|) by adding -0.25 dex and to those with the "diet Salpeter" IMF 



2003h by adding -0.15 dex. 



Even without detailed modelin g, the rest-frame iT-band luminosity for a galaxy is a 



good tracer of its stellar mass (e.g.. 



Bell et al. 



20031). The IRAC 3.6 /im band measures 



rest-frame i^'-band for a galaxy at 0.4 < z < 1.2. Figure \\M shows the i^'-band mass/light 
ratio for the galaxies in four redshift bins. The mass/light ratio is not constant but depends 
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on color (as also found by lBell et al.ll2003l ) and also slightly depends on redshift. Because 
the galaxies are bimodal in color, M/L also shows a bimodality: M/L ~ 1.25 for EC 
galaxies, and M/L ~ 3 for RS galaxies. The different mass-light ratios imply different 
stellar mass limits in the 3.6 /im-selected sample for the two galaxy classes. 

Figure [T5] shows galaxy color versus mass. In general, the most massive galaxies at 
any redshift are the red, quiescent ones. Those are followed by dusty star-forming galaxies, 
which are more massive than the remaining star-forming galaxies in the EC Only the 
dusty star-forming galaxies are massive enough to be direct progenitors of massive red 
quiescent galaxies. Figure US] also shows that there is deficit of lower-mass RS galaxies at 
1.0 < 2; < 1.2. This is not a selection effect: galaxies with M ^ 3 x 10^ Mq would be readily 
detectable regardless of color. The deficit of low mass galaxies on the RS at high redshift 
shows that the most massive RS galaxies formed first, many of them at z > 1.2. The high 
redshift Herschel/ SPIKE sources, representing LIRGs and ULIRGs, tend to be massive, 
indicating that intensive star formation occurs in massive galaxies at these redshifts. 

Figure US] shows galaxy color versus redshift in six stellar mass bins. Most galaxies in 
the highest mass bin are quiescent, indicating that they were already in place at z = 1.2. 
There are a few star forming galaxies in this mass bin as well. On the other hand, the 
most intensive star formation was occurring in massive galaxies with M^, > 10^° M© at 
0.4 < z < 1.2. Figure [T71 shows total infrared luminosities against stellar mass. The plot 
shows that most of the SPIRE sources are LIRGs or ULIRGs, and most are in the greeii 



valley ( iRigopoulou et al. 



Noeske et al. 



(120071 ) 



2012[ ) . Figure [T7] is very similar to Fig. 1 of 
for the visible-selected spectroscopic sample in the same field, showing that more massive 
galaxies have higher total infrared luminosities. Yet the Herschel SPIRE selected galaxies 
are biased towards the high Ljr, comparing with the mean SF R-stellar-mass relat ion in the 



same redshift range for a more sensitive 24/im selected sample ( lElbaz et al 



2003). All of 
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the ULIRGs are more massive than 10^°'^ M©. The mean stellar mass for LIRGs/ULIRGs 
in the sample decreases from lO-'^^'^^ Mq at z = 1.1 to IO-'^'^'^^Mq at z = 0.5, suggesting that 
more massive galaxies migrated to the RS earlier through intensive star formation. 

Figure [16] also shows that quiescent galaxies at z > 1.1 are absent in the 
10 < log(M,/M0) < 10.5 bin and so are those at ^ > 1 in the the 9.5 < log(M,/M0) < 10 
bin. This suggests that quiescent galaxies with low stellar mass did not form until low 
redshifts. Figure [18] shows the RS mass fraction as a function of redshift in each stellar bin. 
Though the mass fraction of quiescent galaxies in 11 < log(M=K/Mo) < 11.5 bin has large 
error bars, it is ~75% and roughly constant from z = to z = 1. This suggests that massive 
galaxies had largely stopped forming stars hj z = 1.2. At 10.5 < log(M*/MQ) < 11, the 
quiesc ent fraction incr eases from 40% at z = 1.1 to 60% at z = 0.6, close to the local value 



(58%, 



Bell et al. 



2003]) at 2 ~ 0. In the 10 < log(M*/M0) < 10.5 bin, the mass fraction 
increases from 10% at 2; = 1.1 to 30% at 2; = 0.5, also close to the local value of 41%. The 
x-intercept of this distribution is at 2; ~ 1.3. Quiescent galaxies in 9.5 < log(M,,/M0) < 10 
are a small fraction of the total at all epochs. The first such galaxies appeared around 
z ^ 1.1, and their numbers have increased to about 25% at 2; ~ 0. 

There are two possible causes for the deficit of red quiescent galaxies in lower mass 
bins in Figure [T8j (1) galaxies with lower mass migrated to the RS later; or (2) quienscent 
galaxies were undergoing either dry or mixed me rging with high ra te to move them quickly 



from lower mass end to high mass end in the RS (lFaber et al 



20071 ). The merging scenario 
requires such a high merging rate in low mass bins th at galaxies just m igrating to the RS 



were immediately undergoing dry/mixed merging. Yet iLin et al.l (120081 ) found that fractions 
of dry and mixed mergers at z~l.l are only 8% and 24% for galaxies with a typical stellar 
mass of 2x10^0 Mq. It is unlikely that dry and mixed merging rate for galaxies with lower 



stellar masses would increase and be the major cause for the deficit of red quiescent galaxies 
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in lower mass at z~l. 

In general, massive quiescent galaxies were already in place hj z = 1.2, and quiescent 
galaxies with 9.5 < log(M^,/Mo) < 10.5 started to form in l.l<z<1.3. This indicates that 
the quenching of star formation, regardless of types, occurs earlier in massive galaxies later 
in galaxies with lower stellar masses. 



4.3. fC-band Luminosity Function and Stellar Mass Functions 

The rest i^-band luminosity function is an important measurement of galaxy population 
evolution. It is equivalent to a stellar mass function but with a better defined completeness 
limit and much less model dependence. This is particularly true for our 3.6 fim selected 
sample. Because the IRAC 3.6 fim probes roughly the rest-frame K-hand for galaxies at 
0.4 < z < 1.2, the 3.6 fim-to-K K-correction is effectively independent of galaxy color, 
permitting an unbiased comparison of the luminosity functions at different redshifts. 
We have derived the i^-band luminosity functions for quiescent, star- forming, and all 
galaxies in the EGS field at 2; = 0.5, 0.7, 0.9, and 1.1 with two meth ods: parametric an d 



non-parametric. The non-parametric method is the famil iar f|Fe 



Bales et al. 



19931 ) ■ and the parametric method is "STY" f lSandage et al. 



ten et al.l 



1976 



19751), which fits 



a Schechter function to the sample. The STY method itself determines only a and 



in the Schechter function. 



variance density estimator (jPavis fc Huchra 



he normalization fact or (p* was derived using the minimum 



19821 ). Both luminosity functions are plotted 



in Figures [20I and numerical values are given in Table 2. 



The major uncertainty in a luminosity function estimate based o n a sample with sna all 



coverage or a single field is cosmic va riance due to galaxy clustering (iHuang et al. 



Somerville 



2004 



Newman et al 



1997 



20021 ) . Our i^'-band luminosity functions were derived from 
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a single l°xlO' field. The maximum angular scale for this sample is 1°, and thus the cosmic 
variance val ues, a,,, at z = .5-1.1 are 28%-36% assuming a general correlation length 
vq = 5Mpc ( ISomervilld 



20041 ). NIR-selected 



that of visible-selected samples: 



g alaxie s hav e a correlation leng th similar to 



Roche et al. 



(120031 ) and 



Oliver et al. 



( I2OO4J ) measured 



correlatio n lengths of 



Recently 



Grazian et al. 



and 6.4 Mpc for K- ari c 



(]2006[ ) and 



Quadri et al. 



IRAC-selected samples, respectively. 



20071 ) measured two-point correlation 



functions for distant red galaxies (DRGs) at z ~ 2 and obtained tq = 7-11 Mpc. The DRGs 
have a typical stellar of 10^^ Mq. If we adopted the largest rg = 13.5, the cosmic variance 
would be 2.5 times larger. However, this large correlation length is not applicable to the 
bulk of our sample but only to the most massive galaxies . The cosmic vari ance can also 



20031 ). The volumes 



be calculated using dark matter variance and bias factor ( iRoche et al.l 
for our sample in 4 redshift bins are ~1.2-2.8 x lO^Mpc^, and within these volumes the 
dark matter fluctuation a£)A4 ~ 0.25. The number densities for L* galaxies in our sample 
are 10^ Mpc~ ^ (Figure l20|) , and therefore the bias for this population at z ~ 1 is 6 ~ 1.5. 



According to 



Roche et al. 



(120031 ) ■ the cosmic variance = &o"dm, and cr^ ~ 0.3-0.4 for the 
galaxies in our sample, consistent with values above. Thus all methods suggest that the 
variance for L* galaxies in this sample is about 30%. 



S-band luminosity functions for galaxie s at 0.4 < z < 1-2 a r e well estab 



in the EGS and a wide range of other flelds (jWillmer et al 



2006 



Faber et al 



ishec 



both 



20071 ). The 



comparison can be used to correct the EGS LF to the global average. The 3.6 /im sample 
is incomplete at faint B magnitudes (Figure fT9l). The 3.6/im limiting magnitude is 23.2, 
corresponding to absolute magnitudes Mk ~ —20 to —21 at 0.4 < z < 1.2. The rest-frame 
B — K colors for galaxies are in the range oi —1 < B — K < 2^ giving a i?-band completeness 
limit of —20.5 to —21.0. The mean density ratio between the EGS and global 5-band LFs, 
taking into account only galaxies brighter than the completeness limit in each redshift bin, 
is {(pegs/4>) = 0.83, 1.7, 1.66, and 1.20 for bins at z = 0.5, 0.7, 0.9, and 1.1 respectively. The 
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standard deviation of these values is 0.38, consistent with the theoretical cosmic-variance 
estimation for this sample. The four values of (0cgs/0) have been applied to the i^-band 
LF and stellar mass function calculated for the 3.6 fim sample. The final LFs are plotted 
in Figure [201 ^ind their values are given in Table 2. These corrected values should be 
representative of the entire Universe, not the EGS field itself. 



The VIMOS VLT Deep Survey (VVDS) flLe Fevre et al.ll2005h provides independent 
i^'-band luminosity functions at 2; = 0.5, 0.7, 0.9, and 1.1. The VVDS itself measured 
redshifts in a wide area also cove r ed by the Spitz e r Wid e- Area Infrared Extragalactic Survey 



(SWIRE) flLonsdale et al 



20041 ). 



Arnouts et al. 



(I2OO7I ) used these data to derive K-hand 



luminosity functions in the same redshi ft bins we have used. Our K-hand LFs are generally 



consistent with those of 



Arnouts et al. 



J, but the VVDS LFs at z = 0.9 and 1.1 are 2 times 
higher than ours at luminosities below L*. We show later that the VVDS sample has higher 
mass densities for star- forming galaxies at 0.4 < z < 1.2 than those obtained in any other 
survey and argue that the excess star-forming galaxies are the cause of their high i^'-band 
luminosity functions at the faint end. Our K-hand luminosity functions a nd the VVDS 



ones with L > L* are higher then the local i^-band luminosity function (IHuang et al. 



2003 



Jones et al. 



20061 ). This is not due to number evolution for massive galaxies because 
massive galaxies can only increase in number towards lower redshifts. Thus it implies an 
luminosity evolution in the i^-band from z ~ 1 to the present epoch. 

Figures EH [22I and [23] respectively show stellar mass functions for all, quiescent, and 
star-forming galaxies. The mass-to-light ratios for quiescent and star-forming galaxies differ 
as shown in Figure [HI in the sense that the mass-to-light ratio for quiescent galaxies is 
roughly three times larger (Figure [T^ than for star-forming galaxies. The survey mass 
limit for quiescent galaxies is therefore correspondingly larger for the fixed 3.6 fim limiting 
magnitude. At the low mass end, star-forming galaxies are the dominant population in 
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the 3.6 /iin sample. Thus the mass hmit for the whole sample should be very close to that 



for star-forming galaxies. Our stell ar mass functions 



consistent with p revious estimates (jProry et al. 



Ilbert et al. 



2003 



br the 3.6 /iin se 



Borch et al. 



2006 



ected sample are 



Bundy et al. 



2006 



20101) but extend below 10^-^ Mq, the lowest stellar mass range so far explored 



at these redshifts. Our stella r mass functions for star-forming galaxies are consistent with 
those in the COSMOS field Jllbert et al.lLoiOl ) except in the bin of OA < z < 0.6, where 
the COSMOS stellar mass function is lower at log(M,,/MQ) < 10.5. The mass functions for 
quiescent galaxies in the EGS and COSMOS fields are consistent with each other at the 
massive end, but our mass functions for quiescent galaxies are higher at the lower mass 
end. Our EGS sample is probably more complete for faint red galaxies because the IRAC 
exposure time is nine times longer than that of COSMOS. 



In comparison with the local stellar mass function (IBell et al.l l2003t iBaldry et al. 



2OI2I ). our stellar mass function shows number evolution since z = 1.1. We also measured 
the stellar mass densities for galaxies with log(M*/M0) > 9.5 in each redshift bin 
(Figure [24)). The total stellar mass density increased by a factor of 3-5 from z = 1.1 to 



0.1 ( Kochanek et al. 



2001 



Bell et al. 



2003 



Driver et al. 



20071 ). In con trast, our mass 



function for 



Driver et al. 



ligh- z star-forming galaxies is similar to local mass functions ( iBell et al. 



2003 



20071 ) (Figure 1231) . implying very weak evolution. The number of star- forming 



galaxies at the lower mass end has increased only slightly since z = 1. The stellar mass 
function for quiescent galaxies shows two different evolutions in Figure [22) the absolute 
number of massive quiescent galaxies and the slope of the low mass end both increase at 
low redshifts. This indicates that massive quiescent galaxies were still forming at z < 1, 
and most low-mass quiescent galaxies formed much later than their massive counterparts. 



The composite stellar mass functions show number evolution since z = 1.1 (Figure [2T]) . 
indicating continuing mass assembly since that epoch. The star-forming galaxy population 
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shows little evolution, and thus the quiescent population is mainly responsible for the 
evolution in the composite mass functions. The quiescent galaxies came from the 
star-forming galaxy population, and the amount of mass assembled in star-forming galaxies 



Ilbert et al. 



also be seen in tht 


i stelL 


ir mass density as a func 


tion of redshift. Several groups ( 


Bell 


et al. 


2003; 


Borch et al. 


2006; 


Franceschini et al. 


2006; 


Abraham et al. 


2007; 


Arnouts et al. 


2007; 



20101) have used stellar mass densities as a function of redshift to demonstrate 



how stellar masses were assembled in galaxies over time (Figure |2^ . Our stellar mass 
densities for both populations are consistent with the other measurements, and indeed 
all surveys yield roughly similar mass densities in the redshift range of 0.4 < 2; < 1.1 
(Figure ^M . Th e mass densities for star-forming galaxies have large scatte r: the VVDS 



has t he highest (lArnouts et al. 



20071 ). and the NEWFIRM has the lowest flllbert et al. 



20101). All surveys yield similar quiescent galaxy mass densities within a factor of 2 in this 



redshift range. This is simply due to better photometric redshifts for quiescent galaxies 
than star-forming galaxies. The overall mass densities increase by a factor of ~4 from 
z = 1.2 to z = 0, and the mass densities for quiescent galaxies also increased by a factor of 
5-8 since z = 1.2. 



5. SUMMARY 

A multi-wavelength study of a 3.6/im selected galaxy sample at 0.4 < 2; < 1.2 with 
/s.e > 2 /iJy in the Extended Groth Strip region has allowed a new suite of photometric 
redshifts based on a neural network method. The distribution of -n — - is Gaussian with 
a ~ 0.025. Pair statistics validates the photometric redshifts even for galaxies too faint to 
obtain spectroscopic redshifts and shows that the redshifts have scatter similar to that for 
the brighter galaxies. The rate of catastrophic redshift failures is about 3.5%. 
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Basic properties measured for the galaxy sample include rest-frame U — B colors, B- and 
X-band absolute magnitudes, and stellar masses for galaxies in the sample. In the redshift 
range available, the IRAC 3.6 /^m band probes the rest-frame near infrared, permitting 
sample selection based on stellar mass and only weakly biased with respect to quiescent 
and star-forming galaxies. Our sample is complete for galaxies with log(M*/M0) > 9.5 at 
0.4 < z < 1.2. We divided the sample into quiescent and star-forming galaxies according to 
their rest-frame U — B color and /24/ /s.e flux ratio. We then derived i^-band luminosity 
functions and stellar mass functions for quiescent, star- forming, and all galaxies in the 
sample z — 0.5, 0.7, 0.9, and 1.1. 

Our conclusions are: 

1. The U — BysB color- magnitude diagram cannot by itself efficiently separate quiescent 
and star-forming galaxies. There are MIPS 24 /xm sources in the red sequence in the 
CM diagram. About half of these sources have f2i/ h.Q > 1, indicating that they are 
dusty star- forming galaxies. The remaining objects have f24/ fs.e < 1, which indicates 
quiescent galaxies with residual star formation or AGN activities. 

2. The U — B-mass diagram shows an absence of quiescent galaxies with log(M*/M0) < 
10 at 1.0 < z < 1.2. The quiescent galaxy mass fractions at low masses 

(9.5 < log(M*/M0) < 10.5) increase with cosmic time from 0%-10% at z = 1.1 
to 25%-35% at z — 0.5. These quiescent galaxies probably first started to appear 
at z ~ 1.3. In comparison, for galaxies an order of magnitude more massive 
(10.5 < log(M*/M0) < 11), the quiescent fraction was already 50% of the total stellar 
mass at z = 1.1, and this fraction has increased only to ~60% at ^ ~ 0. 

3. Throughout the redshift range 0.4 < z < 1.2, ~70% of the massive galaxies 
(log(M*/M0) > 11) are quiescent. These massive quiescent galaxies can have formed 
in one of two ways: either they come directly from star-forming galaxies with similar 
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masses but much shorter hfe times, or they resulted from lower-mass red sequence 
galaxies merging. Some massive galaxies are in a ULIRG phase, giving support to the 
first scenario. 

4. We derived the K-hand luminosity functions and stellar mass functions in four 
redshift bins. The cosmic variance in EGS has been corrected by comparing the 
EGS B-hand luminosity functions with the global B-hand luminosity functions. Our 
X-band luminosity functions suggest there may be a luminosity evolution in the 
K-hand, which needs to be confirmed with larger samples. There is number evolution 
in the stellar mass function, which increases from z — 1.1 to z — 0. Yet the stellar 
mass function for star-forming galaxies shows very weak evolution compared to strong 
evolution in quiescent galaxies. This implies that the amount of mass assembled in 
star-forming galaxies is roughly the same as the amount of mass in galaxies migrating 
from the blue cloud to the red sequence. The faint end slope for the quiescent galaxy 
mass function has increased since z = 1.1, indicating that lower mass quiescent 
galaxies formed much later. The mean mass density increased roughly a factor of 3 
from z — 1.1 to z — 0. 

This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. Support for this work was provided by NASA through an 
award issued by JPL/Caltech. A part of work was done during J.-S. Huang's visit to 
Shanghai Normal University supported by the Chinese National Nature Science foundation 
Nos. 10878003 and 10833005. 
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Fig. 1. — Histogram of number of detections in the multi- wavelength data set. There are a 
total of 18 photometric bands available for the sample. 
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Fig. 2. — IRAC [3.6] — [4.5] colors as a function of redshift. All galaxies in the preliminary 
sample, i.e., before the color cuts of relations 1-3, are shown. The fine dots denote the 



DEEP2 galaxies, and bigger dots denote B]V 


/BX galaxies and L 


(Steidel et al. 


2003, 


2004; 


Shaplev et al. 


2005; 


Huane: et al. 


2005; 



and MIPS 24 /im selected ULIRGs at z ~ 2 flHuang et al 



3Gs detected by 



Rigopoulou et al. 



RAG 



20061 ) 



20091). The line at [3.6] - [4.5] = 



separates galaxies at 2; ^ 1.5: galaxies with smaller redshift have bluer [3.6] — [4.5] color. 
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Fig. 3. — Visible-IR color-color diagram. Points denote all objects selected in the 3.6 /xm 
preliminary sample; red points denote galaxies with spectroscopic redshifts between 1.20 
and 1.30. Lines show the color cuts of relations 2 (lower line) and 3 (upper line). The final 
sample contains only objects in the middle zone of this diagram. 
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Fig. 4. — Visible-IR color-color diagram for objects with redshifts and spectral classifications. 
The locus of the preliminary sample is indicated in gray. Red and blue points indicate galaxies 
in respectively the red sequence and blue cloud in the redshift bins indicated in each panel. 
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Fig. 5. — Montage of visible-IR color-color diagrams as sample selection criteria are succes- 
sively applied. The upper left panel shows all sources in the preliminary sample. The panel 
to its right shows removal of stars via relation 3. The lower left panel shows removal of most 
z > 1.25 galaxies by relation 1, and the final panel shows removal of the rest by relation 2. 
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Fig. 6. — The color-color diagrams for the whole 3.6/im selected sample, the spectroscopic, 
faint (R>24.1) and bright (R<24.1) subsamples. The limiting magnitude for the DEEP2+3 
spectroscopic survey is R=24.1. This diagram shows that faint objects with R>24.1 in the 
sample covers the same color range as the objects with spectroscopic redshifts, and we are 
able to estimate photometric redshifts for these optically faint objects with the spectroscopic 
subsample as the training set in ANNz. 
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Fig. 7. — A direction comparison between photometric and spectroscopic redshifts. Panels 
left to right show the (-Zphot — -2spcc)/(l + ^spec) distributions as a function of redshift for the 
training, validation, and independent validation samples respectively. Upper row shows the 
first ANN2; estimation run, and the lower row the second run with flagged galaxies removed, 
which shows much less systematic deviations. 
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Fig. 8. — Histograms of photometric redshift errors 5zp = (-Zphot — -2spcc)/(l + -^spec)- Panels 
left to right show errors for the training, validation, and independent validation samples 
respectively. Upper row shows the first ANN2; estimation run, and the lower row the second 
run with flagged galaxies removed. Gaussian uncertainty cXp and mean (-Zphot— -^spec) / (l+-2spec) 
for each set are indicated in each panel. 
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Fig. 9. — Histogram of normalized redshift offsets for galaxy pairs selected from the spec- 
troscopic subsample. The red line shows the distribution expected for random, unrelated 
galaxies as calculated by a Monte-Carlo simulation. The insert shows the residual real pair 
distribution after the random distribution of unrelated pairs is subtracted. 
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Fig. 10. — Histogram of normalized redshift offsets for galaxy pairs with one photometric 
and one spectroscopic redshift. The redshift offset distribution for the random, unrelated 
galaxy pair is already subtracted with the Monte-Carlo simulation in this diagram. The red 
dashed line shows a Gaussian fit to the real pair distribution. 
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Fig. 11. — Histogram of normalized redshift offsets for galaxy pairs with one photometric 
and one spectroscopic redshift with the photometric component having R > 24.1. The 
redshift offset distribution for the random, unrelated galaxy pair is already subtracted with 
the Monte-Carlo simulation in this diagram. The red dashed line shows a Gaussian fit to 
the real pair distribution. 
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Fig. 12. — U — B vs Mb color- magnitude diagram for the 3.6 fim selected sample in redsliift 
bins as indicated in the figure. The left panels are for non-dusty galaxies (/24//3.6 < 1); 
and right panels are for dust y galaxies (f^.d/ f?. (\ > !) • The red dots denote Herschel/ SPIRE 



sources selected at 250 /im (IRigopoulou et al. 



20121 ). Most of these sources are LIRGs or 



ULIRGs, and even some of the LIRGs and ULIRGs are in the red sequence. 
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Fig. 13. — Infrared- visible color-color diagram for all galaxies in the 3.6 /xm sample detected 
at 24 [im.. Redshift bins are indicated. The abscissa is the U — B color of each galaxy minus 
{U — B) p as given by equation HI Red dots denote Herschel/SPIRE sources detected at 



250 fim (IRigopoulou et al. 



20121). 
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Fig. 14. — i^-band mass/light ratio versus U — B color. Points denote all galaxies in 
3.6 /xm-selected sample. Redshift bins are indicated in each panel. 
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Fig. 15. — Color versus stellar mass for the 3.6 /xm sample. Redsliift bins are indicated. 
The left panels are for non-dusty galaxies (/24//3.6 < 1), and right panels are for dusty 
galaxies {f'?Al f?..e, > !)• The red dots denote Herschel/ SPIRE sources selected at 250 /xm 



(IRigopoulou et al. 



20121). 
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Fig. 16. — Color versus redshift for the 3.6 /im sample. Stellar mass bins are indicated in 
each panel. The thick lines in the three lower panels indicate redshift limits to detect galaxies 
in these mass bins with tvninal cnlnrs. The limit (z > 1 .51 is outside the diagram in the 
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Fig. 17. — LiR-stellar-mass relation for the Herschel/S PIRE sources in th e 3.6 sample. 



Points denote individual galaxies with Lir values from 



Rigopoulou et al. 



torn . The red 



lines are t he mean SFR-ste llar-mass relation (converted to Lir) for GOODS galaxies at 



0.8<z<1.2 flElbaz et al. 



20071 ). The Herschel galaxies are strongly biased towards the higher 



luminosity side. Redshift bins are indicated in each panel. Histograms show the stellar mass 
distribution of the SPIRE sources, and their median stellar mass is noted in each redshift 
bin. 
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Fig. 18. — Fraction of quiescent galaxies as a function of redshift. Mass range is indicated 
for each panel. Black points denote the average mass fractio n in small r edshif t ranges for 



the 3.6 /xm sample, and red points denote the local value from 



Bell et al. 



(120031). 
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Fig. 19. — i?-band luminosity functions (upper panels) and color-magnitude diagrams (lower 
panels). Redshift ranges are indicated. Points in upper panels denote the lumino sity function 



derive d for the 3.6 fim sample, and lines show the luminosity function derived by 



FaberetaL 



( 120071 ) from a visible sample. Points in lower panels show color and magnitude for 3.6 /im- 
selected galaxies. The absence of points in the lower right of the lower panels shows the 
incompleteness of the 3.6 /im sample for faint, blue galaxies. 
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Fig. 20. — i^'-band luminosity functions. Blue points and solid li nes show re s ults 
this study, while red po ints and dashe d lines show comparisons from 
the leftmost panel and 



Jones et al. 



Arnouts et al 



fl2006h for 



rom 



(120071 ) for the other three panels. Redshift bins are 
indicated in each panel. The lines represent Schechter functions, parameters of which are 
given in Table-El 
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Fig. 21. — Stellar mass functions. Filled circles show results for the 3.6 /xm-selected sample. 
Other symbols show results for comparison samples as indicat e d in t he figure legend. Red 
solid lines indicate the local stellar mass function in 



Bell et al 



(120031). and black solid lines 



Baldry et al 



(120121 ). Vertical 



with arrow are the most updated local stellar mass function in 
dashed lines indicate the completeness limits for the sample in each redshift bin, blue for 
star-forming galaxies and red for quiescent galaxies. Redshift bins are indicated in each 
panel. 
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Fig. 22. — Stellar mass functions for quiescent galaxies. Filled circles show results for the 
3.6 /im-selected sample. Vertical dashed lines show the completeness limit of the sample. 
Open circles and open diamonds show results for quiescent galaxies in the COSMOS field 



( lllbert et al. 



20101) • Redshift bins are indicate d in each p a nel. R ed solid lines show the local 



stellar mass function for quiescent galaxiesin 



Bell et al. 



arrow are the most updated local stellar mass function in 



( 20031). and b l ack so lid lines with 



Baldry et al 



torn . 
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Fig. 23. — Stellar mass functions for star- forming galaxies. Filled circles show results for 
the 3.6 /im-selected sample. Vertical dashed lines show the completeness limit of the sample. 
Open circles and open diamonds s how results for star-forming and starburst galaxies in the 
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Summary of EGS Optical and Infrared Photometry used for Photo-z Estimation 
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Table 2 

EGS X-BAND Luminosity Function 
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Table 3 
EGS Stellar Mass Function 
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EGS Stellar Mass Function for Star-Forming Galaxies 
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